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Abstract 

A novel approach is developed using experimental data 
from the structural testing of a physical systemua identify a. - 
reduced order model and its error for a robust controller de- 
sign. There are three steps involved in the approach. First, an 
approximately balanced model is identified using the Eigen- 
system Realization Algorithm which is an identification algo- 
rithm. Second, the model error is calculated and described 
in frequency domain in terms of the Hoo norm. Third, a 
pole placement technique in combination with a Hoo control 
method is applied to design a controller for the considered 
system. A set of experimental data from an existing setup, 
namely the Mini-Mast system, is used to illustrate and verify 
the approach developed in this paper. 

Introduction 

Constructing a suitable mathematical model for a dy- 
namic system is a major task in the fields of structure analysis 
and control design for vibration suppression. The difficulty 
of estimating and describing the model error,* which is a 
necessity for a robust controller design , 2-12 is even more 
severe than that of constructing a mathematical model. In 
the past, system identification techniques such as the Eigen-... 
system Realization Algorithm (ERA ) 13-18 have been de- 
veloped, using experimental data from the modal testing of 
large space structures, for constructing a mathematical model 
either for a controller design or for modal parameter iden- 
tification including frequencies, damping and mode shapes. 
Yet, the model error issue has rarely been addressed. 

The multi-input/multi-output time domain ERA tech- 
nique was originally developed for modal parameter identifi- 
cation and successfully applied to the modal testing of large 
space structures. The ERA technique 13,1 * uses the singular 
values of a finite Hankel matrix formed by impulse response 
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functions to determine the system order and further reduce 
the noise effect on the realized system model by truncating 
some small singular values. The retained singular values 
are dominated by the system signals and the order of the 
system is determined by the number of the retained singular 
values. On the other hand, the truncated singular values 
are dominated by noise and hence may be used to estimate 
a noise-related model error bound. Defining the model 
error bound to be the maximum truncated singular value is 
equivalent to defining a Hoo error bound for the realized 
system model. Once the Hoo error bound is obtained, a 
Hoo control method 3 may be readily applied to derive a 
robust controller design for the system represented by the 
realized model General speaking, the Hoo control methods 
are optimal frequency domain algorithms which compute the 
feedback control law by minimizing the maximum singular 
value of the system transfer function. It is known that 
the Hoo control methods address the full range of stability 
margin, sensitivity, and robust optimization. This paper is 
motivated by the strong connection between the ERA system 
identification method and the Hoo control methods. It seems 
natural to integrate together the system identification method 
with the control methods so that a robust controller design 

can be achieved directly using the experimental data. 

The objective of this paper is to develop a method using 
experimental data from measurements of a structural testing 
to identify a reduced order model and its error for a robust 
controller design to suppress the vibrational motion of the 
structure. The approach is to integrate the Eigensystem Re- 
alization Algorithm with the Hoo robust control methods. 
This paper is written to outline the fundamental concepts 
from different disciplines and then integrate them together, 
while including advanced materials. In the first section, Han- 
kel singular values of a linear system in frequency domain are 
iji criiccf/l and the relation with the ERA Finite-Hankel singu- 
lar values is described. In the second section, the ERA algo- 
rithm is briefly introduced and the convergence of its model 
realization to the balanced realization 19,20 is described and 
d is cussed. The balanced realization has been widely used for 
model reduction 20 and robust controller designs. In the third 
section, the Hoo control methods are briefly described. In 
the fourth section, the ERA/Z/qo control algorithm is shown 
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including two numerical examples. The ERA reduced real- 
ized model is chosen as the Hqo control design model and 
the estimated model error is used to compute the weighting 
matrices for a robust controller design. In the fifth section, 
an example is given using an existing experimental setup, the 
NASA Langley Mini-Mast system, to illustrate the ERA ///00 
control algorithm. 


where A* is the ith eigenvalue. Now, consider the discrete- 
time state-space model 

z(t+i) = Az(t) + Pm (7) 

K(<) = Cx(t). (8) 

If this system is stable then the discrete controllability and 
observability grammians are defined as 


Hankel and Finite-Hankel Singular Values 
of a Linear System 


/>= £ A k BB*A* k 

CKJkCoo 


Computation of the model error bound is a key element 
in developing robust control systems. The model error bound 
can be expressed and described in many different ways 
including the Hankel singular values of the system transfet— 
function matrix. In this section, the Hankel and Finite* 
Hankel singular values of a linear, time-invariant dynamical 
system are discussed. Consider the linear system 


and 

Q= £ A* k C*CA k 

0<k<oo 

which satisfy the discrete Lyapunov equations 


P - APA* = BB* 


and 

z(t) = Ax(t) -(- B]i(t) (1) 


Q- A*PA = C*C. 


y{t) = Cx{t) (2) 

where x is an n x 1 state vector, u is an m x 1 input vector, 
■g is a / x 1 output vector and A, B, and C are n X n, n x m, 
and /xn system matrices, respectively. The transfer function 
of this system is 


G(s) = C(s/ - A)~ l B. (3) 


The z transform of this discrete system is given as 

G(z) = C(zI-A)~ 1 B 

or in series expansion 

G(z) = fS-z-'- 1 

»=o 


(9) 


( 10 ) 

( 11 ) 


If the system is stable, i.e. all the eigenvalues of A are in the 
left half-plane, then the controllability and the observability 
grammians are defined as 

P= j°° exp(At)BB*exp(A*t)dt 

and 

Q = exp(A*t)C*Cexp(Al)dt (4) 

where P and Q satisfy the Lyapunov equations 
AP + PA* + BB* = 0 
and 

A*Q + QA + C*C = 0. (5) 

Here the superscript * is used herein to indicate the complex 
conjugate of the variable so marked. 

Hankel Singular Values 

The Hankel singular values of this system are defined as 

a i (<7(3)) = (\(PQ)) 1 / 2 , \>Ai + i (6) 


where Y t = CA'B, i = 0, 1, 2,... are impulse response 
functions and usually referred to as the Markov parameters 
associated with the model (A, 6, (7). The Hankel singular 
values of G(z) are defined as 

* i (G(z)) = (X i (PQ)) 1 ' 2 . (12) 

Because a discrete-time model is related to a continuous-time 
model by an invariant step time response, the grammians in 
tire discrete time domain are equivalent to the grammians 
in the continuous time domain. 1 ’ 20 Therefore the Hankel 
singular values for the continuous-time and discrete-time 
models are identical. 

Next it will be shown how the Hankel singular values 
are related to the singular values of a Hankel matrix in 
the discrete-time case. Consider a discrete-time infinite 
Hankel matrix whose (i j)th block is the Markov parameter 
= CA i+ l~ 2 B. The Hankel matrix is written as 



(3 

n 
y * 

... Y-x 
... Y t 


H = W 0 W e = 

iL 

Y r 

• # • 

... Yrl.- J 



i 

* 
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w 0 = [Cr,A , C\...,A* k C\...y 

and , . 

W e = [B,AB,...,A k B,...]. (14) 

It can be shown that P = W e Wg and that Q = W 0 
satisfy Eq. (10). The singular values of H are 

a { (H) = (M/TJ7)) 1 / 2 = (WKWWoWc)) 1 !'. (15) 

In this equation, A i(W*WgW 0 W e ) means that 

W^WlWoWcUi = \Mi (16) 

where u,- is the eigenvector corresponding to the eigenvalue 
Pre-multiplying Eq. (16) by W e and defining x,- = W e Ui 
yields 

W e Vr c W:W 0 Zi = kZi- (17) 

Because Eqs. (16) and (17) produce the same eigenvalues, 
Eq. (15) can be written as 

a^H) = (MVWIOT)) 1 / 2 = (A i(PQ)) l/ *. (18) 

Hence the Hankel singular values of G(z), Eq. (11), which 
is formed from impulse response functions, are the singular 
values of H, Eq. (13). Again note that the Hankel matrix 
H in (13) is infinite dimensional. It is, of course, impracti- 
cal to solve for the Hankel singular values from an infinite 
dimensional Hankel matrix H. If system matrices A, B, 
and C are exactly known, the Hankel singular values can be 
obtained by solving Eq. (12). In practice, however, flexi- 
ble structures are subjected to uncertainties to certain extent 
Experiments are conducted to verify the matrices A, B and 
C. In general, experimental data of inputs and outputs are 
stored in terms of the system transfer function (frequency 
domain) which is equivalent to the impulse response func- 
tions (Markov parameter) (discrete-time domain). Therefore, 
from a practical point of view, the Hankel singular values 
can be approximately computed from a Hankel matrix that is 
formed by a finite number of Markov parameters. The ques- 
tion arises how much error one expects to have if the Hankel 
singular values are computed from a finite-dimensional Han- 
kel matrix. This question will be answered in the following 
section. 


Finite-Hankel Singular Values 

In this section, the singular values obtained from a Finite- 
Hankel matrix are compared with the Hankel singular values. 
First a Finite-Hankel matrix is defined as 


#( 0 ) = 


y 0 y 1 ... 

vi y 2 ... 


,Vr-l Yr 


Y t - 1 

y. 


y-M-2 


(19) 


where Y{ is the Markov parameter. An extended matrix of 
the Finite-Hankel matrix H{ 0) is defined as 


/ Y 0 

tt 

... y,_i o . 




y 3 

... y, o 


II 

Yr-l 

Yr 

• • ■ 
... y r -M- 2 o 



0 

0 

0 0 



; 

; 

• 

*••/ 

It is obvious that 





*(#.) = *(J5T(0)). (2D 


The Hankel matrix H in Eq. (13) can be written as 


H=H?, + K. ( 22 > 


where 


/O ... 0 Y. ...\ 

o ... o y,+i ... 



0 

y+,-1 



(23) 


V 




the left upper r x s blocks of are 0 and the other 
blocks are the same as H. Next it will be shown that 
the Finite-Hankel singular values converges to the Hankel 
singular values. Consider the system with the single-input 
and single-output transfer function G(s) that is real rational, 
scalar-valued, strictly proper and analytic in Re s > 0. Then 
the Markov parameter YJ can be expressed as 


Y { = £c,e**> At (24) 

i=( 

where k is the system order, A ; is the j — th eigenvalue with 
negative real part and A t is the measurement time increment. 
The norm of Yi satisfies the following inequality 

\Yi\ < ce~ xti , U = iAf, i > th (25) 

where c and A are positive real constants, and is a positive 
integer. Using singular value inequalities 21 yields 

am<^ B °r t ) + ^{H‘ rt ). (26) 
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Let r > n,* + 1 and 3 > n,- -f 1 then 

<ri(K,) < W^rtW = 

= f)Yt\ < £ \Yi\ < J2 l ce " Af, 'l (27> 

i=i i=n,+l »=n,*+l 



response functions, 

Y k = CA k B. (29) 

Note that there are many other ways to obtain the impulse 
response functions experimentally. For example, inverting 
a transfer function from the frequency domain to the time 
domain will generate the impulse response functions. The 
algorithm begins by forming the r x 3 block matrix (gener- 
alized Finite-Hankel matrix) 


where e is a positive real number. For any small e > 0, 
there exists an 74 which makes Eq. (27) true. Combination 
of Eqs. (21), (26) and (27) yields 


m = 


Y k Y k+l 
Y k+ 1 Y k+ 2 


Y k+t - 1 
Y k+ . 


<r i (IT)-<Ti{H’ r ,)<<7 i (H?,)<a i (H) + a l (-n‘ r ,) \Y k+r -i Y k + r ... n+r+,-2 


(30) 


or 

Ci(H) - e < Ci(H( 0)) < + e. (28) 

Equation (28) shows that the Finite-Hankel singular values 
converge to the Hankel singular values. The proof can be 
easily extended for the multi-input/multi-output system. 

The finite-Hankel matrix is the basis for the Eigensystem 
Realization Algorithm (ERA) which can be used to identify 
a system model for modal parameter (system eigenvalues) 
identification or controller designs. The Markov parameters 
Yi are the impulse response functions which may be obtained 
experimentally. 

Eigensystem Realization Algorithm (ERA) 
and Model Reduction 

The problem of constructing a suitable mathematical 
model for a dynamic system is a major task to a control 
engineer as well as a structural engineer. For structural anal- 
yses and controller designs experimental data from structural 
testing are used to either verify an analytical model or di- 
rectly determine a mathematical model The ERA is an algo- 
rithm which computes a mathematical model directly from 
experimental data. In this section, the basic ERA formula- 
tions are briefly discussed and the relation between the ERA 
realization and the balanced realization is derived. Note that 
balanced realization has been frequently used to model re- 
duction. 


The Finite-Hankel matrix H(k) can be expressed as 

H(k) = V r A k W, (31) 

where 

j, W, = [B AM ... A- x &\, 

V r and W, are the observability and contro liability matrices, 
respectively. Assume that there exists a matrix n , such that 

W.tfVr = I n (32) 

where I n is an identity matrix of order n. Equations (31) 
and (32) imply 

HWlfHiO) = VrW'rfVrW, = V T W, = H( 0). (33) 

Thus it is the pseudo inverse of H( 0). A solution to H ‘ 
can be obtained as follows. Use singular value decomposi- 
tion to the matrix H( 0) 

#(0) = UZV* (34) 


C 

CA 

\CA r ~ 


Basic ERA Formulations 

The ERA algorithm uses the singular value decomposi- 
tion of a Hankel matrix to determine a system model under 
test Consider the discrete-time state-space model as given 
by Eq. (8). Let the input ji be the impulse input at an ini- 
tial time such that its components satisfy u,(0) = 1, (i = 
l,...,m) and %(k) = 0, (k = 1,...). The time domain de- 
scription is given by the Markov parameters 14 , i.e. impulse 


where the columns of U and V are orthonormal and E is 
diagonal with positive elements [<Ji , o 2 , ...,<7r,]. The pseudo 
inverse can be computed as 

H 1 = [V][E _1 (7 t ]. (35) 

Define 0 m as the null matrix of order m, J m an identity 
matrix, ££ = From Eqs. (31), (32), and 
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(35), a minimal realization can be obtained from 

Y[k + l) = EqH[k)Ep = Ejv r X k W a E p 

= E^V r W a H V r X k W,H V T W,Ep 

= EjUE l l 2 [Z- l l 2 U T H(l)V'L- X l' 1 \ k T>l 2 V T Ep. 

q (36) 

The triple \z- x I 2 u t h(i)vz~ x I 2 , z x / 2 v t e p , eJu e 1 / 2 ] is a 
minimum realization. This is the basic ERA formulation. 

For a finite dimensional and linear time-invariant system, 
an exact system realization can be obtained by the ERA algo- 
rithm from noise-free measurements. 14 The ERA algorithm 
is accurate and efficient, particularly for low noise levels, 
and produces a minimum order realization. It is a power- 
ful identification algorithm. However, if significant noises 
are present in the measurements, caution must be taken to 
identify a proper order for the system model. 

Relation Between ERA and Balanced Realizations 

In this section, the relation between ERA and Balanced 
Realizations is derived. Examination of Eqs. (31) and (34) 
jnriiratpn that the observability matrix and controllability 
matrix can be written as 


written as 

Ct = wjw 0 

= yJvt + (A r ) T v r T VrA r + (A 2r ) T v r T v r A 2r + ... 

= z + (A r ) T rA r + (A 2r ) T EA 2r + ... 

and 

p = W c wJ 

= w t wj + \'w 3 w? (A*) T + A 2a w,wJ (A 2, ) T + ... (40) 
= z + A s e(A*) t + A 2 i E(A 2s ) T + ... 

If the system is stable, all the eigenvalues of A should be 
inside the unit circle, Le., W(A)\ < 1. Therefore, for a 
stable system, there exists a value of n such that |\(A n )| is 
less than z for a given z. Equation (28) implies 

lim E = Eoo (41) 

r— *oo, 8—>oo 

where Eoo is the E in Eq. (34) when r and s approach oo, 
which, in turn from Eq. (40), yields 

lim P = lim E = Eoo 

r— *00, i— ►<» r— *00, t—*oo 


Vr = tfE 1 / 2 and W, = E^V 7 . 
The eq uali ty iPU — I — V*V implies that 
V? Vr = Y}I 2 U t U'£}I 2 = E 


and 

W,Wj = e^v^Ve 1 / 2 = E. 
Combination of Eqs. (14) and (3 1) yields 


/ 


c 

CA 


\ 


Wo = 


CA '- 1 

CA' 


CA 2 '- 1 

CA 2 ' 

i / 



lim 0 — lim E — Eoo • (42) 

r— *oo, #— >oo r— ►oo, *— *oo 

Thus the ERA realization converges to the balanced realiza- 
tion. So far the realization order is assumed to be identical 
to the true one. If the realization system order 7 is lower 
than the true one n, Eq. (36) implies that C^ is the / x 7 
left block of C„. A, is the 7 X 7 left-upper block of A n 
(38) and B-, is the 7 X m upper block of B n . Since the ERA 
realization [A», B n , C n ] converges to the balanced realiza- 
tion [Anfls.Ci]. A,. B n and G, should converge to the 
corresponding blocks of A., £4 and C h , respectively. From 
Eq. (42) and the balanced realization 19 ' 20 , one obtains 


lim P 7 = Ej, 11m — EJ, (43) 

♦00 , j— > oo 1 r~*oo, t—+oo 


(39) where and Q 1 are the controUability and observability 
grammians of [A,, B^, C^] and Eoo is the 7 X 7 left-upper 
block matrix of Eoo- The 7-order ERA realization converges 
to the 7-order balanced realization. It is concluded that 
the reduced ERA model converges to the reduced balanced 
model. 


and 

W. = {B AB...A-'B A‘B...A’-'B A U B..) 

= [W. A‘W, A U W....]. 

The observability and controllability grammians can thus be 


Hoo Robust Controller Designs 

The Hoo control methods are optimal frequency domain 
algorithms 2-4 , which address the full range of stability 
margin, sensitivity, and robust criteria. In this section, a 
brief description of H°° control methods is given. 
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A block diagram of the Hoo design algorithm is presented 
in Fig. 1. The closed-loop transfer function matrices (ui to 
three outputs yj, u 2 and y) are expressed as 


S(s) = (I+L(s))- 1 

(44) 

R{s) = F{s)(I+L(s))- 1 

(45) 

T(s) = Z(s)(/+I(s))- 1 

(46) 


where L(s) = G(s)F(s), G(s) is defined in Eq. (3), and 
F(s) is the controller transfer function matrix. The matrices 
S(s) and T(s) are known as the sensitivity and complemen- 
tary sensitivity, respectively. The singular values of these 
matrices can be used to quantify the stability margins and 
performance of the system. In the Hoo controller design, 
the weighting function matrix W 2 (jw) is used to weight the 
sensitivity so that 

<i(S(jw)) < \W[\jw)\ (47) 

where d( ) means the maximum singular values of ( ). The 
quantity |W / J _ 1 (ju>)| is the desired disturbance attenuation. 

The plant uncertainty can be considered as additive 
uncertainty Ax or multiplicative uncertainty Ajvf. The 
stability margin of the additive uncertainty must satisfy 

o(A A (jw))o(R(jw)) < 1. (48) 

If the weighting function matrix W 2 (jw) is chosen such that 

\W 2 (jw)\ > o(A A (jw)) 

and 

\W 2 {jw)\d(R{jw)) < 1 (49) 

then the inequality Eq. (48) holds. The disturbance attenua- 
tion and additive stability margin designs can be combined 
into a single Hoo problem as 


IPai**i ll°° — 1 

(50) 

where 

T - ( W ' S ) 

lyiUi ~ J 

| (51) 

and this combination is called the Hoo 

additive perturbation 

control design. Similarly, the Hoo multiplicative perturba- 

tion control design can be derived as 


IPviUllloo < 1 

(52) 

where 

T - f WlS ) 

*vi»i - I w 3 T ) 

(53) 


A detailed examination of Hoo design methods can be found 
in Ref. 7. In the next section, a combined ERAJHoq 


control algorithm is developed and discussed including some 
numerical examples. 

ERA/ /Too Control Algorithm 

This section shows a control design algorithm which 
combines the features of both ERA and Hoo control methods. 
The flowchart of this algorithm is shown in Fig. 2. 

To use the ERA/Hgo algorithm, a set of experimental 
input and output data is required. The model error between 
the reduced ERA model and the full ERA model is used for a 
Hoo controller design. The full ERA model is defined as the 
model of maximum order which can be obtained by the ERA. 
For noise-free measurements, if not impossible in practice, 
the full ERA model is the exact system model. For noisy 
meas ur ements; the order of the full ERA model is equal to 
either the number of columns or the number of rows in the 
Finite-Hankel matrix whichever is less. The reduced ERA 
model is the model reduced from the full model such as the 
triple [A,, B-j, C n ] shown in the last paragraph of the ERA 
Section. 

It has been shown that the ERA models converge to 
the balanced models. The 7-th order reduced balanced 
model G~j(s) of an n-th full order system G(s) satisfies the 
following error bound 1 

a(G(jw)-G k (jw))<2 £ <7, (54) 

where ) is the maximum singular value of ( ) and <7; is 
the ith Hankel singular value of G(s). This equation pro- 
vides the Hoo error bound for the reduced model. There 
are two different Hoo control designs, namely the additive 
perturbation control design 11 and the multiplicative pertur- 
bation control design. 11 Numerical examples are used in the 
following to illustrate the ERA/ Hoo control algorithm for the 
two different designs. 

ERA/ Hoo Additive Perturbation Control Design 

In this design the closed-loop transfer function must 
satisfy the inequality 

||T„„|| < 1 (55) 



and u is the input, y is the output, and W\ and W 2 are 
appropriate weighting functions. The sensitivity matrix S 
tends to be dominated by the sensitivity of low frequency 
modes (or dominant modes). An appropriate weighting 
function W 2 is the inverse of the desired sensitivity. The 
weighting function matrix W 2 is chosen to penalize the 
control action at high frequencies (or residual modes) as well 
as incorporate other robustness constraints. One of the design 


6 


goals is to assure that unmodelled high frequency dynamics 
will not destabilize the closed-loop systems. Some examples 
of choosing weighting functions can be found in Refs. 10 
and 11. 

Let the system impulse response measurement be gener- 
ated by the following equation 

4 

y(U) = ^OiexpiQt^sin^iViU) (57) 

t=i 

where the system parameters are listed in Table 1. 


Figure 8 shows the Bode plots of the weighting matrices. 
The Bode plot of the closed-loop transfer function T yii is 
shown in Figure 9, in which the gain margin is 6.0 db and 
the phase margin is 54.2 degree. 

The ERA/TToo control algorithm provides robust control 
designs direcdy from experimental input and output data. 
First, the ERA algorithm is used to compute an appropriate 
reduced balanced model and its error in terms of the /Zoo 
norm. Then, the model error norm is used to select the /Zoo 
control weighting matrices for a robust controller design. In 
the next section, the ERA/Zfoo algorithm is applied to an 
experimental setup, i.e. the Mini- Mast system 22 , for further 
illustration of the concept developed in this paper. 


TABLE 1 



i 1 

2 

3 

4 

a,* 1.0 

0.2 

0.04 

0.01 

Q -0.2 

-0.3 

-0.4 

-0.5 

Wi 3.0 

5.0 

7.0 

8.0 

For illustration, the measurements are 

assumed 


noise free. The Finite-Hankel matrix size is chosen as 
40 x 45. Figure 3 shows the singular-value Bode plots 
of the full model transfer function G, the reduced model 
transfer function <3 7 and the model error A G -G - G 7 . 
Here G n of order 4 is the transfer function obtained by the 
ERA algorithm. The model error A G is an additive model 
error. The additive robustness weighting function and 
the sensitivity weighting function Wi are chosen as 

0.005s + 1 0.04s -I- 0.005 

Wl ~ 0.4s + 0.01’ 2 0.0013 + 1 

Figure 4 shows the Bode plots of A G and W 2 , whereas 
Fig. 5 shows the Bode plots of W\ and Wj. The Bode plot 
of Tyu is shown in Fig. 6, in which the gain margin is 5.2 
db and the phase margin is 56.6 degree. 


ERAJHoo Multiplicative Perturbation Control Design 

In this design, the closed-loop transfer function must 
satisfy 

HTy.ll < 1 (58) 



and Wi and W s are appropriate weighting functions. The 
system measurements and the Finite-Hankel matrix are the 
same as in the above example. Figure 7 shows the design 
robustness W s G n and model error. The design weighting 
matri ces are chosen as 


_ 0,005* -t- X 
Wl “ 0 . 2 * + 0 . 01 ' 


w 2 


= 1.0 x to -7 , W3 = 


0.04* + 0.01 
0.001* + 1 ' 


ERAJHoo Application:* The Mini-Mast System 

In this section, the ERAJH^ algorithm is used to design 
robust controllers for the Mini-Mast system 22 with two in- 
puts (torgue wheels) and two outputs (displacement sensors). 
Mini-Mast is a 20-meter-long, deployable/retractable truss 
located in the Structural Dynamics Research Laboratory at 
NASA Langley Research Center. Mini-Mast was designed 
and built to high standards typical of spacecraft hardware 
which was constructed from graphite-epoxy tubes, titanium 
joints, and precision fabrication techniques. 23 It is used as 
a ground test article for research in the areas of structural 
analysis, system identification, 24 and control of large space 
structures. The sinusoid inputs with the frequencies close 
to the system natural frequencies are used to excite the sys- 
tem. The system has five modes less than 7 Hz and the 
lowest frequency is about 0.86 Hz. The design strategy is as 
follows: 

1. A time domain least squares technique 23 is used to 
compute Markov parameters from input and output 
measurement data. 

2. A Finite-Hankel matrix is formed from the Markov 
parameters and the ERA is used to compute a reduced 
model and its model error. 

3. A Hoo control method is applied to design a robust 
controller. When applying the ERA, the discrete 
time increment should be adjusted to compute an 
accurate model within the frequency range of interest 
Figure 10 shows the singular value Bode plot of the 
ERA model transfer function. The singular value 
Bode plots of the model error (solid lines) and the 
rfr rignfrd additive robustness (dashed lines) are shown 
in Fig. 11. 

In MATLAB 7 , the Hoo control computations are 
performed by using the recent Glover-Doyle “2-Riccati” 
formulae. 12 After the Hoo control design is executed, the 
closed-loop system keeps the open-loop stable complex 
eigenvalue pairs and reflects the open-loop unstable complex 
eigenvalue pairs. Because the identified open-loop model of 
the Mini -Mast system has small positive damping (system 
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inodes) or negative damping (noise modes) for most complex 
eigenvalues, the Hoo closed-loop system will then have low 
damping. To improve this condition, a multivariable bilinear 
transform 8 can be used to shift the open-loop eigenvalues 
before applying the Hoo control methods. The process of 
the bilinear transform is described below: 

1. The system model and the design weighting functions 
are shifted with a real value a. 

2. The Hoo design methods are used to design a con- 
troller for the shifted system. 

3. Either the closed-loop system is shifted back by —a 
or the controller is used in the original nonshifted 
system. 

4. Check the singular values of the closed-loop transfer 
function matrix T„ u . Eq. (56), which .should -be. less 
than one. 

In the process of applying the bilinear transform, several 
problems listed below are found. 

1. The original system has a Hoo control solution 
with appropriate design weighting functions, but the 
shifted system with the shifted weighting functions 
does not always have a Hoo control design solution. 

2. The singular values of the closed-loop system in- 
crease significantly when the imaginary axis has a 
positive shift 

To improve the above problems, a different technique which 
is equivalent to a pole placement technique is developed 
yielding freedom for assigning each closed-loop eigenvalue 
independently. Let the system state matrix be diagonalized 
such that 

A = 'fU'Ir ( 60 ) 

where $ is the eigenvector matrix. Let the closed-loop state 
matrix be represented as 

A = *(A+ Al)*" 1 (61) 

where Al is a semipositive diagonal matrix. The design 
strategy is as follows. 

1. Use A as the state matrix for the Hoo controller 
design. 

2. Apply the controller to the original system. 

3. Check whether this control design satisfies the Hoo 
constraint equation, Eq. (55). 

Figure 12 shows the singular value Bode plot of a Hoo 
design transfer function T yu with a specified matrix Al 
which is not shown. The matrix Al provides more than 
10% damping for the first mode. Next, the computational 
simulation is used to check the performance of the above 
design. The solid lines in Fig. 13 are the first 5-sec 
experimental outputs. In these plots, represents the 
ith output with the jth input The outputs from the ERA 
reduced model using the same inputs as the experimental 


data are plotted in dashed lines. The inputs are turned off 
after 5 seconds. The open-loop and closed-loop responses 
are compared in Fig. 14. The solid lines in Fig. 14 are the 
open-loop free decay responses, whereas the closed-loop free 
responses are plotted in dashed lines. The first closed-loop 
mode has more than 10% damping. 

In this ERAJHoo control design, the following items are 
observed. 

1. This algorithm can be used to design a robust con- 
troller directly from a set of experimental input and 
output data. 

2. The controller order is the sum of the model order 
and weighting function order. Economical controllers 
may be obtainable through the use of appropriate 
model reduction techniques. 

3. If the identified low frequency modes are somewhat 
accurate, the robustness weighting can cover signifi- 
cant high frequency model errors. 

4. The additive matrix Al can be used to assign the 
closed-loop dampings and frequencies for the modes 
of interest. 

5. An analytical model is not required. 

Concluding Remarks 

In this paper, the Eigensystem Realization Algorithm is 
integrated with die Hoo control methods to develop a tech- 
nique for controller designs of flexible structures. The tech- 
nique starts with a set of input and output data obtained from 
vibration tests of a physical system. The set of data is then 
used to identify a reduced model and assess its error bound. 
The reduced model describes the physical system, whereas 
the error bound characterizes the system and measurement 
uncertainties. Both reduced model and error bound are used 
to design a robust controller to attenuate the vibrational mo- 
tion of the physical system. Since experimental data are 
used, actuator dynamics, sensor dynamics and filter dynam- 
ics are included in the identification of the reduced model 
and thus in the controller designs. This approach is believed 
to provide the reduced model with fidelity. Numerical sim- 
ulations for control of an experimental setup (the Mini-Mast 
system) indicate that this approach is very promising for use 
in control of flexible structures. Note that experimental mod- 
els are used in these simulations. The controllers which are 
successfully demonstrated in the simulations will be tested 
in the Mini-Mast system to verify the approach developed in 
this paper. 
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Figure 1 Weighted sensitivity and complementary sensitivity problem 



Figure 2 ERA/H^ Design Flowchart 
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Figure 4 Design Robustness and Model Error 
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Figure 6 Closed- loop Transfer F\mction T yu 
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Figure 9 Closed-loop Transfer Function T yu 
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Figure 10 System Model 
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